1 



a^ 
a^ 

> 
o 

o 

(N 



> 

o 



Exotic Stochastic Processes from Complex Quantum Environments 

Dimitri KUSNEZOV^ *, Aurel BULGAC'^t^ and Giu DO DANG^ 

'^Center for Theoretical Physics, Sloane Physics Laboratory 
Yale University, New Haven, Connecticut 06520-8120, USA 

''Department of Physics, University of Washington 
P.O. Box 351560, Seattle, Washington 98195-1560, USA 

'^Laboratoire de Physique Theorique, Universite de Paris-Sud 
Batiment 211, 91405, Orsay, FRANCE 



Abstract 

Stochastic processes are shown to emerge from the time evolution of complex quantum systems. Using parametric, 
banded random matrix ensembles to describe a quantum chaotic environment, we show that the dynamical evolu- 
tion of a particle coupled to such environments displays a variety of stochastic behaviors, ranging from turbulent 
diffusion to Levy processes and Brownian motion. Dissipation and diffusion emerge naturally in the stochastic 
interpretation of the dynamics. This approach provides a derivation of a fractional kinetic theory in the classical 
limit and leads to classical Levy dynamics. 
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1. Introduction 

The understanding of how stochastic processes 
emerge from classical dynamical systems is closely 
related to classical chaos. Often one finds that the 
dynamics is non-Gaussian, displaying either en- 
hanced or dispersive behavior ||l|. One can find 
such behavior in the interaction of slow and fast 
degrees of freedom in many-body systems , in 
tracer diffusion in turbulent backgrounds such as 
the atmosphere, or random potentials, and many 
more . The wide variety of processes which ex- 
hibit anomalous behavior in the transport has led 
to a variety of theoretical efforts, including frac- 
tional extensions of kinetic theory |^-^ , random 
walks in random potentials Q, power law noise 
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in generalized Langevin equations stochastic 
webs [|| and Levy walks and flights |||,||. While 
the common thread to these approaches are gen- 
eralizations of Brownian motion known as Levy 
stable laws (discussed below) , there is no common 
theoretical foundation. For instance, fractional 
kinetic equations are postulated in such a man- 
ner as to provide the desired diffusion through 
scaling arguments. 

In turning to the quantum theory, we might ask 
whether it is possible to realize anomalous diffu- 
sion or Levy stable laws in the time evolution of 
quantum Hamiltonians ||l^. If classical chaos is 
the origin of the stochastic processes in classical 
dynamical systems, it is natural to ask whether 
quantum chaos can be the source of quantum 
stochastic processes. Starting with the quantum 
counterpart of classical chaos, namely random 
matrix theory (RMT), we consider whether the 
time evolution of such Hamiltonians can generate 
stochastic processes. We will see that not only is 
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this possible, but depending on certain properties 
of the quantum Hamiltonian, one can reahze a full 
range of stochastic processes, including those of 
Levy. Further, fractional kinetic theory develops 
naturally in the semi-classical limit. 

To understand diffusion observed in many sys- 
tems which have {R^{t)) ^ f , where 7 7^ 1, 
generalizations of Brownian motion have been 
sought. Levy processes provide such an approach. 
In studies of the extensions of the central limit 
theorem, P. Levy found a continuous class of non- 
Gaussian processes that satisfy the same funda- 
mental equation that gives rise to the theory of 
Gaussian processes ||T^. In one dimension, the 
Levy stable laws have the form 



dk 
2^ 



exp{ikx- A\k\"} (1) 



where < a < 2 and A oc t. Gaussian processes 
correspond to the case where a = 2. The Levy 
distributions are scale invariant, 



(2) 



where for A — I wc drop the superscript: 
C^{x) = Ca{x). The scale invariance indicates 
that the trajectory followed by the random pro- 
cess will not be dominated by one characteristic 
scale, resulting in a self-similar behavior. With 
the exception of the Gaussian (a = 2), all the 
Levy stable laws have infinite second moments 
For certain values of a, it is possible to com- 
pute the inverse Fourier transform in (1) [jl^ . 

2. Quantum Chaotic Environments 

We would like to understand the diffusion of 
a quantum particle interacting with a quantum 
chaotic environment. Random matrix theory pro- 
vides a convenient description of quantum chaotic 
systems, but as it contains no scales or physics, 
we must develop the notion of a random environ- 
ment. Since our approach is Hamiltonian based, 
we envision the particle interacting with a com- 
plex quantum system, which might have a non- 
trivial density of states, and which is inherently 
chaotic in the sense of RMT fluctuations of the 
matrix elements. We will assume that the Hamil- 
tonian Hf, which describes this 'environment' is 



time-reversal invariant, so that it is a real, sym- 
metric matrix. (This serves only to simplify the 
notation). As the background is chaotic, it is not 
necessarily thermal. 

The Hamiltonian for the chaotic environment 
plus interaction will have the form 



He = ho{x,p) + hi{X,x,p). 



(3) 



where (a;,p) are the coordinates and momenta of 
the environment, and (X, P) are those for the test 
particle, /iq is the Hamiltonian of the environ- 
ment and hi is the interaction with the test par- 
ticle. It is convenient to choose the fixed basis of 
/iQ to describe the matrix elements of hi. In this 
basis, we denote 

/io I n) =e„ I n), (n = l,...,iV), (4) 
so that the matrix elements of H^. are 



[HeU =e.Aj + [hi{X)], 



(5) 



There is now control over the density of states 
of ft-o by suitably choosing the £i to describe the 
system of interest. We would like to introduce 
a parameter related to the average level density 
of /iQ as /3 = 1/T = d£np{e)/de, or equiv- 
alently, p{e) = po ex-p{(3e). We will see that this 
derivative of the level density naturally plays the 
role of temperature in the quantum dynamics in 
certain cases. We now choose hi to be chaotic, 
which allow us to apply the Gaussian orthogonal 
ensemble (GOE) to the matrix elements. In this 
case, the chaotic properties of the interaction be- 
tween the environment and the test particle are 
built into the correlation function (second cumu- 
lant): 



[hi{x)],j[hi{Y)]ki) = g^,{x - y)A 



ijkl • 



(6) 



Here A^-fe; = [5ik6ji + Sa^jk]: and aU other cu- 
mulants vanish. A convenient realization of the 
correlation function is ]T3|-p^: 



2TT^p{£i)p{ej) 

X exp 



2kq 



G 



The essential elements of this function are the 
following. We expect due to selection rules that 
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hi is not a full matrix, but most likely banded 
in some way. Here the band width is described 
through kq. As the energy increases, the level 
density typically increases rapidly, so one expects 
the average interaction matrix elements to re- 
duce, which is accounted for by the fac- 
tors. One also expects He to decorrelate in X on 
some length scale Xq (NB This is not generally 
the scale on which the eigenvalues of He{X) vary 
as suggested by the avoided level crossings). Fi- 
nally, G(x) = Gi-x) = G*{x) < 1, G(0) = 1, and 
the overall strength is known as the spreading 
width. Before we turn to the dynamics, let us 
consider the statistical nature of hi. 

3. Short Distance Correlations 

An interesting consequence of the statistical 
measure leading to (6) is that the correlation 
function given by the two-point function must 
decorrelate no faster than quadratic jl^] . In many 
cases, the essential characteristic of G{x) is gov- 
erned by its short distance behavior. On short 
distance scales we approximate: 



Gix)^l-Ca\x\" + 



(8) 



For systems with smooth parameter dependence, 
one has a = 2. A Dyson process on the other 
hand has a ~ 1 [p^ . Similarly, disordered sys- 
tems are often modelled with a — 1, for instance 
when one uses {V{x)V{x')) ~ exp[— |a; — x'\]. 

From the requirement that the probability dis- 
tribution of matrix elements be always 
a positive definite function and using the theo- 
rems of Bochner, it can be shown that the corre- 
lator G{x) has to satisfy the following restriction 
at short distances 



< a < 2. 



(9) 



Hence, we obtain a finite range of values a which 
are allowed for parametric random matrix ensem- 
bles. 

As the position X of the slow particle 
changes, the instantaneous energy levels En{X) 
of [hi{X)]ij change. Using the above expresion 
for the correlator G{x) one obtains that the aver- 
age fluctuations are 



The energy-spacing fluctuations have a behavior, 
which is similar to a Levy process characterized 
by the diffusion constant Da- The character of 
these fluctuations in the eigenvalues En{X), in- 
dicated by a, will be seen to be related to Levy 
distributions, which describe the time evolution 
of the density matrix for a particle evolving in 
this chaotic bath. 

4. Quantum Levy Processes 

We now derive the dynamics of a quantum par- 
ticle coupled to this quantum chaotic environ- 
ment. We use the Hamiltonian 



H,j{X,P)^S,. 



p2 

+ UiX) 

2M ^ ' 



+ (11) 



We will be interested in the diffusive and dissipa- 
tive behavior of the test particle induced by the 
environment, so we will neglect the effects of the 
potential U{X). The evolution equation for the 
density matrix of the test particle for this class 
of RMT Hamiltonians has been derived recently 
using influence functional techniques |p^ . It is 
given in a high temperature expansion, and up to 
o(/3) has the form: 



ap(xy,t)^rp£_^ 

dt 1 2M 2M ^ ' ^ ' 



AXqM \ Xo 
X-Y 



G 



X, 







- 1 



{Px - Py) 
p{X,Y,t). 



(12) 



{[E,,{X) - En(Y)f) = D^\X -Y\' 



(10) 



To extract the dynamics analytically, we pass to 
the weak-coupling limit []l^ , in which we keep 
only the leading order term in the correlation 
function. This corresponds physically to the case 
where the test particle does not provide signifi- 
cant feedback to the environment. In this limit, 
we have G{x) = 1 — |a;|" so that G'{x) represents 
— Q!sign(a;)|a;|"~"'". 

In the absence of a potential U{X), and when 
the level density is constant in the region of in- 
terest {(3 — 0), these equations are readily solved 
for the test particle density matrix p{X,Y,t). In 
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the variables r ^ {X + X')/2, s ^ X - X' , the 
solution is 



p(r, s, t) = / dr' I (r', s - — 



X exp 



zA:(r — r') 



-kt/M 



(13) 



The initial density matrix at i = is denoted by 
piX,Y,0)^PoiX,Y). 

The probability of finding the test particle at a 
position X at a certain time t, denoted P{X, t), is 
given by the diagonal elements of the density ma- 
trix: P{X, t) = p(X, X, t) = p{r, s = 0, t). This is 
readily computed from Eq. (20) as: 



2h^ 8Ma^ 



(14) 



We have assumed that the test particle is initially 
in a Gaussian wave-packet 



i^iX) = 
PQ{r,s) = 



na 



2] 1/4 



2TTa 



exp 



--L(4r2 + .^) 



(15) 
(16) 



Comparing to Eq. (1), we see that Eq. (21) is a 
Fourier transform of a Gaussian and a Levy pro- 
cess, which can be expressed as the convolution 
of the two distributions, 

P{X,t)^ J dX' C<'\X')cf\x - X'), (17) 



with 

a{t) = 



{a + l)h \MXo 



bit) 



a 
~2 



8Af2cr2 



(18) 



(19) 



It is interesting to note that these processes are 
Markovian, having no memory effects which are 
responsible for the anomalous character. The dy- 
namics of these processes are now labeled by a. 



which was shown in Section 3 to be in the range 
< a < 2. Hence, Eq. (17) shows us that 
the short distance decorrelations of the adiabatic 
states of ife are directly related to the type of 
stochastic process which develops in the quantum 
time evolution of the test particle. 



5. Turbulent like Diffusion 

For the cases a < 1 and a > 1, the character of 
these Levy processes has been discussed in px[] . 
When a — 2 and /3 = 0, P{X, t) is the convo- 
lution of two Gaussians, and as a consequence is 
Gaussian as well. The test particle in this case 
exhibits turbulent diffusion, with 



{X'^;t) = j dX X^p{X,X,t) 



4M2cr2 3M2X2 



(20) 



(21) 



Here one can see that the initial width at t — Q 
is due to the Gaussian initial condition, and the 
second term reflects the natural spreading of a 
wave-packet in the absence of an environment. 
The character of the dissipative contribution 
which arises from the environment is typical of 
turbulent backgrounds. The range of times un- 
der which we expect such diffusion to hold has 
been discussed in Ref. ||l^. In Table 1 we sum- 
marize selected properties of the RMT influence 
functional and its stochastic behaviors. 



6. Fractional Kinetic Theory 

One of the phenomenological approaches to 
anomalous diffusion is to postulate a Fokker- 
Planck (FP) equation which supports the types 

of diffusion seen in complex classical systems. 
To this end, a number of extensions of the FP 
equation have been proposed in which the spa- 
tial and/or time derivatives are replaced with 
derivatives of fractional order. Fractional deriva- 
tive/integrals are typically defined through the 
use of integral transforms |Q . In one of the many 
realizations, for example, one can define I?" to 
be a derivative (integral) of order a when a > 
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{a < 0) as follows. For a > let 
1 



dy{x-yr-'f{y), (22) 



where a > is real, f{x) is an arbitrary function 
and r(a) is the F-function. For f{x) = x^, where 
^ > is real, this gives 



r(/i + i) 



+a 



(23) 



Clearly for a — n a positive integer, this is just 
the n-th iterated integral of f{x). For general 
a > it is a suitable definition for the fractional 
integral of order a. For negative values of a, let n 
denote the integer part of a = —a. For this range 
of a we take 



In this case 



F(n-a + l) 



(24) 



dy {x-yr-^f{y). 



(25) 



which is what one might expect for a derivative 
of order a. Fractional derivatives can be related 
to anomalous diffusion by noting that a diffusion 
process with (a;^) = Dt"' can be derived from 
a fractional Fokker-Planck equation through the 
use of scaling arguments Q . 

There are many realizations of the fractional 
FP equations to physical systems, including ap- 
plications from turbulence to diffusion in vis- 
coelastic of porous media . Generally they are 
of the form ||,| 



VlP{Q,t) 



V^_Q{A{Q)P{Q,t)) 
+ Vq {B{Q)P{Q,t)) 



(26) 



with fractional derivatives in space and/or time. 
The powers <5, /x, v are chosen to reproduce 
anomalous diffusion through scaling formulas 
such as ~ V , where 7 is a function of 5, fi, v. 
But there are many limitations at the moment to 
this type of phenomenological approach. First, 
one must assume the existence of the coefficients 



A, B and so forth, which are now fractional mo- 
ments. Further, since there is no microscopic un- 
derpinning, extensions to higher dimensions be- 
come tenuous since many more coefficients are 
needed, whose origin is then not understood. 

By starting with the quantum problem of a 
test particle in the chaotic quantum environment, 
we have introduced certain 'microscopic' scales 
which describe the interactions between the sys- 
tems. If we now take a semi-classical limit of our 
influence functional, we will see that its classi- 
cal analog is a fractional Fokker-Planck equation. 
Further, the quantum distribution functions are 
now solutions (in this limit) of that equation. 

The fractional Fokker-Planck equation is de- 
rived by taking the Wigner transform f{Q,P,t) 
of the density matrix p{X, Y, t) 



f{Q,P,t) = 



dR 

2TTh 



exp 



iPR 



(27) 



The classical transport equation is now defined 
for f{Q,P,t) by Wigner transforming the evolu- 
tion equation (19). To leading order in h, 



dfiQ,P,t) 
dt 



dR 



■ exp 



iPR 



R 



xp{Q + 



R 



ijhXoasign{R) 



R 



R 



Xo 



(28) 



We now use the Reisz fractional integro- 
differential operator which is defined as ]2^Jl0[| 



(-Ap)t/(p) = .F-i|xr.F/(p), 



(29) 



where Ap is the Laplacian (here with respect to 
the momentum P), and J-' is a Fourier transform 
from P to X. For our purposes we take Dp = 
(-^/?i)"(-Ap)"/^ since = d^f/dP^. We 

also use Dp = {-i/h)°'J'-hign{X)\X\°'J', which 
has the property l3p[P/] = d{Pf)/dP. The 
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result is a fractional Fokker-Planck equation in 
phase space, which is also a fractional extension 
of Kramers equation: 

df{Q, P,t) P dfiQ, P, t) dU{Q) df{Q, P, t) 
dt M dQ dQ dP 



„{i?P '[P/(Q,P,<)] 
2 (*?i)"i?^[/(Q,P,i)] 



^ 7. 

2TM 



(30) 



ah 



Here T = 1//3 is the temperature and the opera- 
tor and the generalized friction coefficient is 



7q 



2MX 



(31) 







The solutions to this equation follow directly from 
the Wigner transform of our quantum solutions 
(20) and (24). 

7. Conclusions 

The quantum dynamics of test particles in com- 
plex quantum environments was studied using in- 
fluence functional methods. The stochastic char- 
acter of the dynamical evolution follows from the 
specific properties of the chaotic/random matrix 
environment. Markovian dynamics ranging from 
Levy to turbulent and to Brownian diffusion are 
possible. In the classical limit, one can derive a 
fractional Fokker-Planck equation, which reduces 
to the well-known Kramer's equation for the spe- 
cial case of a complex environment. This provides 
a connection between quantum chaos (RMT), 
stochastic processes and fractional kinetic theory, 
which can be used to better understand the ori- 
gins of anomalous diffusion in complex classical 
and quantum systems. 
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Table 1. Selected dynamical limits of the ran- 
dom matrix influence functional of Ref. [[l5| . 



dcLSsic8/l Limit 


a = 2 
a<2 


Kramer's Equation [|l5| 

Fractional Fokker-Planck Equation; Levy Processes 


U{x) = 


/3 = 0, a = 2 


Turbulent Diffusion ||l|l 




(3 = 0, a<2 


Quantum Levy Processes [|l0| 




/5 7^ 0, BM limit 


Quantum Brownian Motion with (x^) ~ 2Dt 




f3 ^ 0, Xo^l 


(x'^) — 2Dt but non-Maxwellian distributions ||l5|]. 


Weak Coupling 


/3 > 0, a = 2 


Caldeira-Leggett Influence Functional jl^,^ 



